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Multiplicity distribution in small bins is studied within the Ginzburg-Landau description for 
first-order quark-hadron phase transition. Direct comparison of the distribution with a Poisson one 
(with the same average) is made. Dynamical factor dq for the distribution and ratio Dq = dq/di 

QQ , are studied, and novel scaling behaviors between Dq are found which can be used to detect the 

^\ • formation of quark-gluon plasma. 
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description 

: 

' It is well-known that high energy heavy-ion collision is the unique way to study the vacuum properties of quantum 
chromodynamics (QCD) in laboratory. In the collisions a hot new matter state, quark-gluon plasma (QGP), might be 
' formed, and the system will cool with its subsequent expanding and will undergo a phase transition from deconfined 
^ , QGP to confined hadrons. Since only the final particles in the collisions are observable in experiments, one may be 
^ ' asked to search for the signals about the phase transition from only those particles. Contrary to the final state leptons, 
(y-^ ' which are produced mainly in the early stage of the new matter state and do not participate the strong interactions, 
T-H ' hadrons are emitted out in the freeze-out process and may carry some relic information about the evolution of the 
' system and the interactions. Up to now many possible signals about the phase transition, such as J/^* suppression, 
. strangeness enhancement, etc., have been proposed theoretically. But all of them can be reproduced from conventional 
physics. Thus, more further studies about the signal are needed. 

In this Letter, we study the final state hadrons for the signal based on the Ginzburg-Landau description of the 
phase transition. Speaking generally, multiplicity distribution of final state hadrons is one of the most important and 
most easily accessible experimental quantities in high energy leptonic and hadronic collisions. The study of the global 
distribution, from the earlier well-known KNO scaling and its violation [1, 2] to a recently novel scaling form [3], can 
O 1' give us a lot about the dynamical features of the multiparticle production processes. Local multiplicity distribution has 
O . also been studied for many years in terms of a variety of phase space variables [4] , and substantial progress has been 
1-^ ' made recently in deriving analytical QCD predictions for those observables [5] . Based on assuming the validity of the 
^ , local parton hadron duality hypothesis, these predictions in [5] for the parton level can be compared to experimental 
k>l data. A global and local study of multiplicity fluctuations [6] shows, however, that those theoretical predictions have 
rN [ significant deviation from experimental data. The significant deviation of theoretical predictions from experimental 
■ data indicates that we know only a little about multiparticle production processes since the hadronization process in 
■ ■ ' soft QCD is far from being understood, and phenomcnological investigation may still be valuable. 

It is well-known that an important feature associated with critical phenomena is large fluctuations. For multiparticle 
production processes the scaled factorial moments [7] is the most effective way to extract dynamical fluctuations. In 
fact, the moments are averages of some quantities over multiplicity distribution. Then one can see that the dynamical 
information for the processes can also studied in terms of the distribution. In particular, one may expect that whether 
the colliding system undergoes a quark-hadron phase transition will be reflected in the distribution. In this Letter 
we try to investigate multiplicity distribution in some small two-dimensional kinematic bin (which can be rapidity 
and transverse momentum or azimuthal angle, for example) in high energy heavy-ion collision processes. Because the 
final state hadrons may carry some relic information about their parent state, the investigation of the multiplicity 
distribution may be interesting and useful for probing the formation of QGP. As is shown in a previous work [8], 
the study can reveal more information than that of factorial moments can since the average over multiplicity in 
the calculation of the moments may conceal some useful information. In this Letter we are limited to discussing 
multiplicity distribution under the assumption of a first-order phase transition, within Ginzburg-Landau description 
for the phase transition. The multiplicity distribution for the case of second-order quark-hadron phase transition 
has been studied in [8], in which two new quantities are introduced to measure the influence on the distribution of 
dynamical fluctuations, and interesting scaling results are shown. It should be pointed out that the scaled factorial 
moments are studied within the same description for phase transitions by many authors for both second-order [9, 10] 
and first-order [11-14] phase transitions, and a universal scaling exponent ~ 1.30 is given in [9,10,13,14]. 
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In Ginzburg-Landau theory, the multiplicity distribution turns out to be a Poisson distribution if the field is pure 
coherent. Conversely, the distribution turns into a negative binomial distribution if the field is completely chaotic. In 
reality, one can assume multiparticle production arising from a mixture of chaotic and coherent field, so the multiplicity 
distribution in real processes is not a Poisson one nor a binomial one, and the deviation of the distribution from a 
Poisson one is due to dynamical fluctuations. The real quantity concerned with the dynamical process is the relative 
deviation of Pg (the probability with q hadrons in a bin) from its Poisson counterpart Pq. Thus studying the relative 
deviation may reveal features of dynamical mechanisms involved. As introduced in [8], the relative deviation of Pq 
from its Poisson counterpart pq can be measured by ratio dq = Pq/Pq- The ratio dq can be called dynamical factor, 
since it is 1.0 unless there are dynamical fluctuations in the process. Dynamical fluctuations is shown to be exist if the 
ratio is far from 1.0, either much larger or much smaller than 1.0, for some q. For the definition of dq to make sense, it 
is necessary to let the mean multiplicity (n) for Pg and Pg be the same which will be given below. In Ginzburg-Landau 
description for first-order phase transition the distribution Pg is given in [9] 

Pqi6) = / VcPpqiS^ I <p He-^I^l , (1) 



where Z = J P^e--^!"*! is the partition function, Pq{s) is the Poisson distribution with average ,s, Pq{s) = s"^ /q\ exp(— s), 
and F[<p] the generalized free energy functional for the first-order phase transition in a two-dimensional bin 6^ 

Fld^] = 52 [a I |2 I |4 +c I </. |6] , 

with & < 0, c > 0, and a cc T — Tc representing the temperature. For more information about the discussion of the 
fimctional, see [8]. For different sets of parameters the system can be in the quark phase or hadron phase, see [11-14] 
for detail. In real experiments the temperature at which hadrons are emitted from the source is unknown and may 
be different from event to event. So we treat a as a free parameter and discuss only fluctuations in hadron phase with 
a < in the following since in the quark phase only a few hadrons can be produced through fluctuations. Prom the 
distribution of Eq. (1) one gets the mean multiplicity for a < 

with Hn{u,v) = dy y'^ exp{—y^ + uv'^y + vy"^) and s = (J^/c)-^/^ representing the bin width 5, B = b/^/c. This 
expression is different from that for second-order phase transition case in that the mean multiplicity in flrst-order case 
depends on three (instead of two) variables, the bin width s, the temperature | a |, and B. The mean multiplicity in a 
bin is proportional to s at small phase space bin widths, thus it can be very small. For very small bins the distribution 
must be concentrated around Pq, and both Pg and Pq for g > 1 must be very small. So a direct comparison between 
Pq and pq could induce large uncertainty. This demands that the bin width in real experimental analysis should be 
large enough to ensure the mean multiplicity in the bin not too small (larger than 0.5, say). Of course, smaller bins 
can be used if a precise determination of both Pg and Pg can be obtained from data with high statistics. Naively, one 
may expect the distribution to approach to a Poisson one in the small s limit. In the Ginzburg-Landau description 
for flrst-order phase transition, it is not the case. In the small s limit, the mean multiplicity is approximately 

sr(2/3)/r(i/3), 

" - q\ r(l/3) ' ^« - g! U(l/3) ) ' ^ ' 

Thus except Pi , which approaches to pi in the small s limit, Pq ^ Pg for q > 1 \a the small s limit. Because of 
the normalization of both Pg and Pg, the dynamical factor dg must be larger than 1.0 for some q and less than 1.0 
for some other q if there exist dynamical fluctuations. Using the Ginzburg-Landau description for first-order phase 
transition one can easily get 



-yfsB) gg(| a I /B\~^sB) 
Ho{\ a I /B2, ^^sB) HU\ a \ /B^, ~^sB) 



To see the dependences of (n) and dg on the bin width s, we graph (n) and Indq as functions of — Ins in Fig. 1 
for B = —1.0, I a |= 1.0 and 2.0, respectively. One can see that for larger s (or smaller — Ins) the mean multiplicity 
is larger and that in dg is larger than 1.0 (Inrf^ larger than 0.0) for some q but less than 1.0 for some other q. For 
smaller s (or larger — In s) simple increasing behaviors can be seen for all dq with the increase of — In s. This means 
that the dynamical influence can be observed more easily in smaller windows in the kinetic regime. This can be 
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understood once one notices the fact that different dynamical fluctuations may offset each other and become less 
obviously observable in large bins. In the limit of s 0.0 all dg approach to certain saturation values. Another 
feature for the dependences of (n) and In dg is that the detailed dependences rely on the temperature \ a\, especially 
obvious in larger bin width s. 

It is instructive to notice that the ratio Pg{{n))/pi{{n)) is an increasing function of (n) for g > 1, although Pq{{n)) 
itself exhibits complicated behavior with the increase of (n). Moreover there exists a scaling law between pg 

Pg{{n)) _2'i-^ fp,{{n))Y-' _ 



Prom this scaling law, one may conjecture that it is more interesting to study the dependence on s of 



Dg^^ = ^ (6) 

di Pq/P 



instead of dg and that one may expect similar scaling behavior between Dg when the resolution is changed. One can 
see from below that it is true. Now we turn to study Dq for first-order quark-hadron phase transition. If there is no 
dynamical reason, Pq — pq, one can see that Dg for all q can have only one value, 1.0, no matter how large or small 
the bin width is. So from the range of values Dq takes, one can evaluate the strength of dynamical fluctuations. Dg 
can be expressed in terms of Hn{u, v) as 

Besides the bin width s, there is in last expression two more parameters B and | a \ which is a measure of how far 

from the critical temperature the hadronization process occurs and is unknown in current experiments. First let us fix 
B = —1.0 and let | a |=1.0 and 2.0, respectively. The dependence of Dq on — Ins is shown in the upper part in Fig. 2. 
Now Dq for all q is an increasing function of — Ins. For small — Ins values of Dq depend strongly on parameter | a |, 
but they approach parameter independent values for large —Ins. The similarity among Dq for different q suggests 
that there may exist a scaling law between Dg and D2. The same data Dg for q > 2 are reshown in the lower part 
in Fig. 2 as functions of D2- Power scaling laws between Dq and D2 arc found for both | a |=1.0 and 2.0. For other 
values of I a I the similar power law dependence is checked to be true. Thus one has 

\nDq = Ag + BqlnD2, (8) 

with Aq and Bq depending on | a | for fixed parameter B — —1.0. The fitted results of Aq and Bq from curves in 
lower part in Fig. 2 are shown in Fig. 3 as functions of ln{q — 1) for | a |=1.0 and 2.0. It is obvious that both In 
and InBq have linear dependence on ln{q — 1) for flxed \ a \. Especially, for studying power law, we investigate Bq 
and find that 



Bq = {q-iy, (9) 

with 7 depending on | a | for fixed B = —1.0. For visualization, the linear fitting curves for InBq vs ln(g— 1) are shown 
also in Fig. 3 for | a |=1.0 and 2.0. The slopes for Aq is about twice those for Bq. The exponent 7 increases with 
increasing \ a \. When | a | is zero, corresponding to the case in which hadrons are produced exactly at the critical 
point, numerical results show that 7 is 0.912 for fixed B = —1.0. With the increase of | a | at fixed B, 7 increases 
quickly. For sufficiently large | a |, when the difference between | a | — 1 and | a | can be neglected, corresponding to 
the case in which hadrons are produced at temperature far below the critical point, one finds that Dq Fq, with Fq 
the scaled factorial moment in the same kinematic bin which is given in [10,11] for first-order phase transition as 

Similar relation between Dq and Fq is also true in the small s limit. In these limiting cases, the scaling of Dq is 
equivalent to that for the scaled factorial moments Fq, and one can get the exponent 7=1.4066 for large — Ins [15] 
or large | a | . This exponent is a little different from those obtained from the study of the factorial moments for both 
first- and second-order phase transitions [9, 10, 13, 14] and of the multiplicity difference correlators for second-order 
phase transition [16]. The dependence of 7 on ) a ] is shown in Fig. 4 for fixed B = —1.0. The dependence of the 
exponent 7 on parameter B is also studied. As an example, we fix parameter | a ]=1.0. It can be seen clearly from 
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Eq. (7) that 7 is equal to 1.0 for parameter B = 0.0 which corresponds to the tricritical case. Since we are discussing 
first-order phase transition, B should be negative. For larger —B the exponent 7 is larger, as shown in Fig. 5. For 
other fixed parameter | a \ similar dependence of 7 on _B can also be seen. 

In conclusion, dynamical factors dq and their ratios Dq are studied to describe features of dynamical fluctuations 
in first-order quark-hadron phase transition. In Ginzburg-Landau description for first-order quark-hadron phase 
transition it is found that Dq obeys a power law, Dq cx I?^', with Bq ^ [q — 1)'^ . Both dq and Dq can be obtained 
quite easily in experimental analysis. First one can get Pq by counting up the number of events with exactly q hadrons 
in the bin. pq is of Poisson type and can be calculated from the experimental mean multiplicity (n) in the bin. Simple 
algebras give dq and Dq. The existence of dynamical fluctuations can be confirmed if dq and Dq can take values very 
different from 1.0. The range of values for Dq is a measure for the strength of dynamical fluctuations. The scaling 
between Dq and D2 is a possible signal for the formation of QGP, because up to now no other dynamical reason is 
known to induce such a scaling. The study of Dq should be carried out in real experimental analysis in the near 
future to see whether QGP has been formed in current high energy heavy-ion collisions. As tools to study dynamical 
fluctuations dq and Dq introduced in this paper may also be interesting in experimental analysis of leptonic and 
hadronic interactions without quark-hadron phase transitions. 

This work was supported in part by the NNSF, the SECF and Hubei NSF in China. 
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Figure Captions 

Fig. 1 Dependences of the mean multiplicity in a bin and the dynamical factor dq (9=1, 2, 3, 4, 5, and 6) on the 
bin width - Ins for fixed B = -1.0, | a | = 1.0 and 2.0. 

Fig. 2 Upper figure: Dependences of Dq on bin width — Ins for | a |= 1.0 and 2.0 and B = —1.0 for g=2, 3, 4, 5, 
6; Lower figure: Scaling behaviors between Dq and D2 for B = —1.0, | a |=1.0 and 2.0. 

Fig. 3 Coefficients for the scaling between Dq and D2, InDg = Aq + BqlnD2, as functions of \n{q — 1) for fixed 
B = —1.0, I a 1 = 1.0 and 2.0. Linear fitting curves are shown for Bq = [q — 1)'^ . 
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Fig. 4 Dependence of exponent 7 on | a | for fixed B = —1.0. For large | a |, 7 is about 1.4. 

Fig. 5 Dependence of exponent 7 on —B for fixed | a |= 1.0 for first-order quark-hadron phase transition {B <0). 
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